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INTRODUCTION 

It  has  been  known  for  several  hundred  years  that  the  earth  possesses  a 
planetary  magnetic  field.  It  has  also  been  known  for  over  200  years  that 
there  Is  a  dally  variation  In  the  surface  strength  of  that  field.  The  earth's 
surface  field  also  exhibits  changes  associated  with  solar  flares.  Even  when 
the  sun  Is  "quiet1',  however,  there  Is  a  variation  In  the  earth's  surface 
magnetic  field.  This  dally  variation  at  mid  latitude  Is  quite  regular  and  Is 
typically  represented  by  "S^",  which  stands  for  solar  quiet  dally 
variation.  Typically  at  mid  latitude  this  dally  variation  ranges  from  20  to 
45  nanoteslas  (nT).  This  Is  compared  with  a  mean  surface  magnetic  field  of 
about  50,000  nT.  Near  the  turn  of  the  century,  the  existence  of  the 
Ionosphere  was  Inferred  from  these  magnetic  variations.  It  was  suggested  that 
the  Sq  pattern  was  caused  by  Ionization  In  the  upper  regions  of  the 
atmosphere  being  moved  by  neutral  particles  through  the  earth's  magnetic 
field.  This  eventually  became  known  as  the  Ionospheric  dynamo  theory. 

When  the  existence  of  the  Ionosphere  was  proven  by  reflecting  radio  waves  off 
of  It,  the  Ionospheric  dynamo  theory  became  firmly  entrenched.  With  the 
discovery  of  the  magnetosphere,  however.  It  was  realized  that  other  currents 
flowing  above  the  earth  should  also  contribute  to  the  surface  magnetic  field 
variations.  In  early  work  on  these  magnetospherlc  currents.  It  was  suggested 
that  Indeed  they  do  contribute  to  the  earth's  surface  magnetic  field,  but 
produce  an  Inconsequential  portion  of  the  observed  pattern.  Since  those 
reports,  much  has  been  learned  concerning  the  nature  of  the  magnetospherlc 
current  systems.  The  purpose  of  this  work  Is  to  quantitatively  reexamine  the 
question  of  non-lonospherlc  currents  and  their  contribution  to  Sq.  (It  was 
expected  that  somewhat  different  results  might  be  obtained  because  of  our 
Increased  understanding  of  the  magnetospherlc  currents.)  It  Is  Important  to 
understand  what  this  contribution  Is  since  It  will  put  an  upper  bound  on  the 
neutral  winds  In  the  upper  atmosphere  required  for  the  dynamo  theory.  It  Is 
also  Important  to  quantitatively  determine  the  magnitude  of  the  magnetospherlc 
currents  at  the  earth's  surface  In  order  to  accurately  determine  the  earth's 
main  magnetic  field.  The  accurate  determination  of  the  main  field  (to  within 
several  nanotesla's)  Is  needed  for  the  study  of  the  secular  variation  In  the 
main  field,  and  for  the  accurate  description  of  magnetic  anomalies  In  the 
earth's  crust.  In  addition,  the  development  of  the  capability  for 
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representing  the  contribution  of  magnetospherlc  currents  In  the  earth's 
surface  magnetic  field  should  also  be  helpful  In  the  study  of  magnetospherlc 
dynamics,  since  ground  based  magnetograms  are  routinely  used  to  qualitatively 
describe  magnetospherlc  events. 

HAGNETOSPHERIC  CONTRIBUTION  TO  SURFACE  FIELD  -  DIRECT  INTEGRATION 
One  of  the  problems  with  early  work  on  the  magnetospherlc  magnetic  field  was 
that  models  representing  It  were  meant  to  be  global  -  (l.e.,  useful  over  the 
entire  magnetosphere)  and  little  care  was  given  to  the  accuracy  of  the  field 
from  the  particular  magnetospherlc  current  at  the  earth's  surface.  Therefore 
It  Is  expected  that  significant  errors  persist  In  the  surface  magnetic  field 
predicted  by  these  early  models.  Therefore  the  Blot-Savart  law  was  used  to 
Integrate  over  the  magnetospherlc  current  systems  to  directly  determine  the 
magnetic  field  at  the  earth's  surface.  This  was  done  Individually  for 
magnetopause,  ring  and  tall  current  systems.  The  results  are  shown  In 
Figure  1.  It  Is  seen  when  the  magnetospherlc  magnetic  field  models  are  used, 
that  the  total  contribution  from  the  magnetospherlc  currents  at  the  earth's 
surface  produced  a  day  to  night  variation  of  about  4  nT.  However,  when  a 
direct  Integration  Is  performed  over  the  current  systems,  this  variation  Is 
actually  about  7  nT.  Note  that  these  numbers  refer  to  the  magnetic  field 
produced  In  the  region  of  the  earth  but  assume  that  a  zero  conductivity  earth 
Is  present,  l.e.,  the  field  Is  given  In  free  space.  The  pattern  for  the 
ring,  magnetopause,  and  tall  currents,  and  their  sum  Is  shown  In  Figures  2,  3 
and  4  for  the  north,  east  and  vertical  components  of  the  total  magnetic 
variation  field  employing  direct  Integration  over  the  currents.  In  each 
figure,  the  dally  variation  Is  shown  at  0  ♦  30°,  and  ♦  60°  magnetic 
latitude  and  the  vertical  tick  marks  are  In  5  nT  Intervals.  Note  that  In 
these  figures  a  free  space  (zero  conductivity  earth)  was  assumed. 

However,  several  problems  remain.  First,  although  there  Is  some  uncertainty 
In  the  representation  of  all  three  major  magnetospherlc  current  systems 
(magnetopause,  ring,  and  tall),  It  Is  the  ring  current  that  causes  the  largest 
qualitative  uncertainties.  As  shown  In  Figure  5,  It  Is  seen  that  different 
portions  of  the  Inner  part  of  the  ring  current  can  be  set  to  zero  and  have 
little  effect  on  the  total  ring  current  contribution  In  those  regions  of  the 
magnetosphere  where  It  Is  most  routinely  observed.  However,  these  changes  In 
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the  total  ring  current  strength  dramatically  Influence  the  value  of  the  ring 
current  In  the  vicinity  of  the  earth.  Another  shortcoming  of  work  to  date  Is 
that  the  Blrkeland,  or  field  aligned  currents,  have  not  been  considered.  It 
Is  clear  that  they  make  an  Important  contribution  to  the  surface  magnetic 
field  near  auroral  latitudes. Finally,  the  results  In  Figures  2  through  4 
assume  a  zero  conductivity  and  zero  permeability  earth. 

THE  QUANTITATIVE  DETERMINATION  OF  EARTH  INDUCTION  CURRENTS  AND  THEIR 
CONTRIBUTION  TO  THE  SURFACE  MAGNETIC  FIELD 

It  has  been  known  through  this  century  that  the  earth's  electromagnetic 
properties  cause  the  magnetic  field  from  external  sources  to  be  significantly 
modified  In  Its  vicinity.  Thus  It  Is  expected  that  the  surface  magnetic  field 
produced  by  the  magnetospherlc  currents  will  be  modified  considerably  by  the 
real  earth.  It  was  therefore  necessary  to  quantitatively  study  the  problem  of 
the  currents  electromagietlcally  Induced  In  the  earth  as  It  rotates  under 
these  magnetospherlc  current  systems.  Although  work  on  Induced  currents  and 
their  associated  magnetic  fields  was  not  a  part  of  the  original  work 
statement,  It  was  felt  that  because  of  the  Impact  of  Induction  currents  on  the 
total  magnetospherlc  Sq  value.  It  was  appropriate  to  quantitatively 
determine  their  contribution.  Typically  the  earth  has  been  considered  as  free 
space  or  as  having  Infinite  conductivity.  A  finite  conductivity  earth  has 
been  used  typically  only  In  Induction  problems  where  earth  currents  are 
Important  (e.g.,  In  the  study  of  crustal  anomalies).  Using  the  magnetospherlc 
magnetic  fields  (as  shown  In  Figures  2-4)  as  the  forcing  function,  the  dally 
variation  In  the  Sq  pattern  was  found  for  various  values  of  conductivity  for 
the  earth.  There  Is  no  question  about  either  the  existence  or  the  magnitude 
of  these  magnetospherlc  currents  nor  their  contribution  to  the  pattern 
and  that  the  absolute  minimum  they  contribute  to  the  day  to  night  amplitude  Is 
12  nT. 

We  have  solved  the  three-dimensional  problem  of  the  magnetic  fields  Induced 
Internal  and  external  to  a  sphere  having  constant  conductivity  which  Is 
rotating  In  the  presence  of  a  static  external  magnetic  field.  The  theory  for 
determining  these  Induced  fields  Is  developed  In  the  Appendix. 
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In  order  to  solve  for  the  Induced  magnetic  fields.  It  Is  necessary  to  express 
the  forcing  external  fields  as  an  Internal  spherical  harmonic  expansion.  This 
has  been  done  and  the  coefficients  of  the  expansion  are  given  In  Table  1.  The 
magnetopause  expansion  Is  taken  from  Mead  (1964),  assuming  a  10  Rg 
magnetopause  standoff  distance.  The  rlr.g  and  tall  expansions  come  from 
magnetic  field  models  developed  at  MOAC  for  zero  ■tilt''  (Olson  and  Pfltzer, 
1974). 

The  Induced  magnetic  field  found  outside  the  sphere  Is  represented  by 
coefficients  of  an  external  spherical  harmonic  expansion.  The  coefficients  on 
Table  2  give  these  coefficients  assuming  the  sphere  has  a  conductivity  of 
.0065  mho/meter.  This  Is  the  "preferred*  surface  conductivity  given  In 
Alldredge  (1977),  and  since  conductivity  Increases  dramatically  with  depth 
below  the  surface,  this  model  probably  underestimates  the  actual  induced 
magnetic  fields.  The  conductivity  dependence  of  the  total  surface  field  Is 
examined  In  Figures  6  through  9.  Figure  9  shows  the  north-south  field 
component  at  the  equator.  Notice  that  the  field  variation  Is  Increased  using 
the  "preferred"  conductivity  by  SOX  over  the  non-conducting  value.  The  same 
holds  true  at  30°  latitude  for  both  the  north-south  component  (Figure  8)  and 
the  east-west  component  (Figure  9).  The  radial  component  variation  (shown  In 
Figure  7)  Is  reduced  to  one-third  of  Its  non-conducting  value  assuming  the 
same  conductivity. 

The  value  .0065  mho  Is  most  commonly  used  to  represent  the  conductivity  of  the 
earth's  crust.  Since  the  earth's  conductivity  Increases  with  depth,  this 
value  somewhat  underestimates  the  ability  of  the  earth  to  Induce  currents  and 
shield  Its  Interior  from  the  magnetospherlc  magnetic  field.  In  absolute 
terms,  the  day-night  variation  of  the  forcing  magnetic  field  Is  7-3/4  nT  at 
the  equator.  With  the  "preferred"  conductivity  of  Alldredge,  an  Induced  field 
Is  produced  to  Increase  this  variation  to  11-3/4  nT,  Implying  that 
magnetospherlc  sources  contribute  a  large  part  of  the  ~  30  nT  day-night 
variation. 

The  Sq  pattern  produced  by  the  3  magnetospherlc  current  systems  and  their 
associated  Induced  currents  Is  shown  for  the  north  east  and  vertical 
components  of  the  surface  field  at  0,  ♦  30°  and  ♦  60°  latitude  In 
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EXTERNAL  MAGNETIC  FIELD'S  INTERNAL  HARMONIC  EXPANSION  FORCING  FUNCTION* 


***Fit  of  zero  tilt  current  wire  model  magnetic  field. 
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Figures  10-12  with  the  earth's  conductivity  set  at  .0065  mho.  It  Is  seen  that 
the  amplitude  of  the  variation  Is  Increased  to  about  12  nT  when  the  earth's 
electrical  conductivity  Is  taken  Into  account  (as  represented  by  the  dashed 
lines).  The  "free  space  earth"  values  are  represented  by  the  solid  lines. 
(Mathematical  considerations  for  the  representation  of  the  Induced  earth 
current  systems  are  given  In  the  Appendix). 

Note  that  as  the  conductivity  goes  from  zero  to  Infinite,  the  phasing  of  the 
variation  pattern  changes  such  that  the  maximum  contribution  occurs  before 
noon  for  all  finite  conductivities,  but  occurs  exactly  at  noon  for  both  zero 
and  Infinite  conductivity.  It  Is  also  worth  noting  that  the  phasing  problem 
has  been  used  to  suggest  the  magnetospherlc  currents  do  not  make  an 
appreciable  contribution  to  Sq.  The  argument  goes  that  because  the  maximum 
In  typically  occurs  after  noon,  and  also  typically  occurs  after 
hemispheric  solstice,  only  lags  In  thermospheric  heating  associated  with  the 
Ionospheric  dynamo  currents  can  explain  the  observations.  However,  phasing 
caused  by  the  earth's  finite  electrical  conductivity  and  also  phasing  In  the 
annual  variation  pattern  caused  by  the  fact  that  the  solar  wind  Is  not 
Incident  upon  the  geomagnetic  field  from  the  optical  solar  direction,  both 
suggest  that  the  magnetospherlc  contribution  possesses  the  appropriate 
patJ  ..  (The  Walters  effect  suggests  that  the  free  stream  solar  wind  makes 
an  angle  of  about  4  degrees  to  the  optical  sun-earth  direction  and  It  has  been 
argued  that  when  the  solar  wind  flows  through  the  bow  shock,  the  angle  between 
Its  flow  direction  and  the  optical  sun-earth  line  may  exceed  10  degrees). 

VARIABILITY  IN  THE  S  PATTERN 

- q - 

In  order  to  quantitatively  study  the  day  to  day  variability  In  S^,  It  was 
necessary  first  to  develop  a  data  base  containing  ground  based  magnetometer 
observations  and  Information  on  the  variability  of  solar  wind  parameters.  Of 
course,  both  data  sets  must  be  for  the  same  period  of  time.  We  were  able  to 
obtain  IMP-B  solar  wind  data  for  1977  and  1976  from  NASA  Goddard.  We  were 
also  able  to  obtain  mean  hourly  values  of  the  three  components  of  the  vector 
magnetic  field  observed  at  about  two  dozen  mid-latitude  magnetic  observatories 
during  the  same  two  year  time  period  from  the  World  Data  Center  In  Boulder, 
Colorado.  We  had  hoped  with  these  data  sets  to  be  able  to  show  that  the  day 
to  day  variability  In  the  earth's  surface  magnetic  field  Is  caused  primarily 


16 


NORTH-SOUTH  FIELD 


(D^seiouDU)  HldOS-HiaON 


17 


12.0  18.0  24.0 

MAGNETIC  LOCAL  TIME 


EAST-WEST  FIELD 


o 


(o^seiouou)  1S3M-ISU3 


18 


MAGNETIC  LOCAL  TIME 


VERTICAL  EIELD 


(D^se^ouou)  iyOU^3A 


9 


0.0  6.0  12.0  18.0  24.0 

MAGNETIC  LOCAL  TIME 


by  changes  In  the  magnetopause  magnetic  field.  We  used  solar  wind  data  to 
describe  changes  In  the  magnetopause  standoff  distance  (and  also  the  strength 
of  the  magnetopause  current  system).  The  change  In  the  magnetic  field 
produced  by  the  magnetopause  currents  at  the  earth's  surface  at  noon  In 
reponse  to  solar  wind  pressure  Is  shown  In  Figure  13.  The  solid  curve  shows 
simply  the  magnetic  field  from  the  magnetopause  currents;  the  dashed  lines 
shows  the  total  magnetic  field  when  earth  currents  are  Included.  It  Is  seen 
that  as  the  magnetopause  standoff  distance  changes  from  about  10-1/2  down  to  7 
earth  radii,  the  contribution  from  the  magnetopause  currents  to  the  earth's 
surface  magnetic  field  changes  by  approximately  60  nT.  Thus  we  hoped  to  find 
this  variation  present  In  the  day  to  day  changes  In  the  surface  magnetic 
field.  The  data  from  the  ground  observatories  was  restricted  to  quiet 
magnetic  conditions,  as  evidenced  by  low  Kp  values.  Further,  In  order  to 
remove  the  contribution  from  the  ring  current,  the  Ost  value  was  removed  from 
the  data  (at  the  earth's  surface,  Dst  directly  represents  the  magnetic  field 
produced  by  the  ring  current).  The  horizontal  component  of  the  field  at  San 
Juan  Is  shown  at  noon  for  quiet  days  with  the  ring  current  contribution 
removed.  The  expected  strong  dependence  on  solar  wind  dynamic  pressure  was 
not  detected,  as  Is  shown  In  Figure  14.  A  similar  plot  Is  shown  for  Newport 
In  Figure  15.  The  dependence  on  the  magnetopause  currents  would  tend  to  make 
the  data  points  lie  along  the  curve,  which  shows  the  computat1ona'lly 
determined  contribution  of  magnetopause  currents  at  the  earth’s  surface  at 
noon  as  a  function  of  dynamic  solar  wind  pressure.  Note  that  the  curve  could 
be  moved  horizontally.  We  also  note  that  If  such  a  correlation  can  be 
filtered  out  of  the  data.  It  could  be  used  to  uniquely  determine  the  baseline 
of  Sq.  We  are  unaware  of  any  other  way  to  do  this  quantitatively,  thus  our 
study  of  the  dependence  of  the  surface  magnetic  field  on  solar  wind  variations 
will  continue  because  an  accurate  determination  of  the  earth's  main  field  Is 
Important  for  many  theoretical  and  commercial  reasons. 

OVERALL  ASSESSMENT  OF  HA6NET0SPHERIC  CONTRIBUTION  TO  THE  SURFACE  MAGNETIC  FIELD 
It  Is  difficult  to  assess  precisely  what  fraction  of  the  observed  Sq  pattern 
these  magnetospherlc  currents  produce.  This  Is  because  there  Is  at  present 
some  Indecision  concerning  the  magnitude  of  the  Sq  currents.  Classically, 
the  mid  latitude  pattern  exhibits  a  peak  to  peak  variation  of  about  25 
gammas.  However,  recent  Investigations  of  Sq  suggest  that  the  pattern  for 
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very  quiet  days  during  the  last  decade  exhibits  a  day  to  night  variation  of 
approximately  45  gammas.  Thus  the  magnetospherlc  contribution  of 
approximately  12  gammas  should  produce  somewhere  between  one-fourth  to  almost 
one-half  the  current  of  the  pattern.  There  Is  also  an  uncertainty  In  the 

value  of  the  magnetopause  ring  and  tall  currents  that  have  been  used.  These 
have  been  average  values  for  each  current  system  as  defined  by  magnetospherlc 
configuration  for  quiet  times  over  long  periods.  However,  It  Is  known  that 
the  magnetosphere's  current  systems  are  quite  variable  as  they  respond  to 
large  fluctuations  In  solar  wind  and  Interplanetary  field  parameters.  It 
would  be  appropriate,  wherever  possible,  to  determine  the  strength  of  the 
magnetospherlc  current  systems  for  the  same  days  that  were  used  to  derive  the 
observed  patterns.  Clearly  the  magnetospherlc  values  presented  here  are 
baseline  values  and  a  detailed  Investigation  day  by  day  (where  solar  wind  and 
Interplanetary  field  data  are  available)  should  yield  somewhat  larger  values 
for  the  magnetospherlc  contribution.  We  therefore  remain  unconvinced  of  the 
Ionosphere's  dominant  role  In  producing  surface  magnetic  field  variations 
during  quiet  magnetic  conditions. 

Since  the  contribution  we  are  now  finding  for  these  currents  Is  approximately 
2-1/2  times  larger  than  reported  earlier.  It  Is  reasonable  to  assume  that  this 
will  also  be  true  for  the  annual  variation  produced  by  the  magnetospherlc 
currents.  This  would  suggest  that  the  magnetospherlc  currents  then  produce 
approximately  one-fourth  the  observed  8-10  nT  seasonal  variation. 

Work  done  over  the  past  two  years  as  part  of  the  Coordinated  Data  Analysis 
Workshops  (CDAWs)  has  shown  that  although  the  nominal  position  for  the 
magnetopause  Is  between  10  and  11  earth  radii  along  the  sun-earth  line.  It  Is 
sometimes  compressed  to  well  within  geosynchronous  orbit.  The  variability  In 
magnetospherlc  features  during  the  July  29,  1977  event  was  extensively  studied 
In  COAW-2.  The  magnetopause  location  ranges  from  about  6  to  over  15 
In  just  over  half  a  day  (see  Figure  16).  Compression  of  the  magnetosphere  Is 
responsible  for  the  well  studied  "sudden  commencement  phase"  of  the 
geomagnetic  storm.  Using  pressure  balance  formalism  It  Is  easy  to  show  that 
an  Inward  motion  of  the  magnetopause  of  1  earth  radius  will  Increase  the 
magnetopause  contribution  to  Sq  by  approximately  35  percent.  Thus  Increases 
In  the  surface  field  for  magnetospherlc  currents  as  large  as  50  and  even  100 
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percent  are  not  uncommon.  Studies  of  the  solar  wind  suggest  that  there  are 
definitely  extended  periods  where  the  solar  wind  pressure  Is  considerably 
enhanced,  and  yet  the  solar  wind  flow  Is  steady.  Thus  It  Is  clearly  possible 
that  the  magnetospherlc  contribution  to  Is  larger  than  reported  here.  It 
Is  also  clear  that  the  magnetospherlc  currents  must  make  a  significant 
contribution  to  the  day  to  day  variability  In  the  pattern.  We  expect 
that  solar  wind  pressure  variations  between  quiet  days  accounts  for  a  large 
part  of  the  range  In  the  strength  of  the  observed  pattern  for  those 
days. 

In  order  to  accurately  determine  the  main  magnetic  field  (as  has  been 
attempted  recently  with  MAGSAT  data).  It  Is  necessary  to  have  a  capability  for 
representing  the  magnetospherlc  currents  In  real  time.  A  problem  with 
attempting  to  average  magnetospherlc  data  over  some  Interval  of  time  (e.g.  an 
epoch  used  for  the  MAGSAT  analysis)  Is  that  the  dependence  of  the  Sq 
contribution  from  the  magnetopause  location  Is  a  nonlinear  function  of 
magnetosphere  size.  Thus  the  use  of  magnetospherlc  models  to  accurately 
determine  the  surface  magnetic  field  requires  the  availability  of  time  varying 
representations  of  each  of  the  magnetospherlc  current  systems.  These  have 
been  developed  and  are  now  being  used  as  part  of  the  CDAW  process.  They  may 
also  be  used  In  conjunction  with  MAGSAT  data  sets  In  an  attempt  to  accurately 
determine  features  of  the  earth's  Internal  magnetic  field. 


SUMMARY 

The  purpose  of  this  work  has  been  to  reexamine  the  role  that  the  magnetosphere 

plays  In  producing  magnetic  variations  at  the  earth's  surface.  Prior  to  the 

discovery  of  the  magnetosphere.  It  had  been  assumed  for  half  a  century  that 

the  magnetic  variations  at  the  earth's  surface  during  quiet  times  (the  Sq 

pattern)  were  caused  by  currents  flowing  In  the  Ionosphere.  Some  of  the  more 

astute  "founding  fathers"  typically  referred  to  the  source  of  S  as  caused 

Q 

only  by  an  equivalent  "overhead  current  system"  Instead  of  specifically 
suggesting  that  all  of  the  contribution  was  caused  by  the  Ionosphere.  Early 
work  on  the  magnetosphere  suggested  that  It  produced  an  Inconsequential 
portion  of  this  Sq  pattern.  However,  the  early  magnetospherlc  models  were 
developed  to  represent  the  magnetic  field  through  the  entire  magnetosphere  and 
were  not  constructed  for  accuracy  at  the  earth's  surface.  Because  of  this  and 
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our  much  better  understanding  of  the  magnetospherlc  currents  fifteen  years 
later,  another  look  has  been  taken  at  the  contribution  of  the  magnetospherlc 
currents  to  Sq.  So  far  It  has  been  found  that  the  magnetospherlc 
contribution  to  S^,  as  represented  In  terms  of  a  day  to  night  range.  Is 
approximately  12  nT  when  the  earth's  finite  conductivity  Is  taken  Into 
account.  (Although  the  earth  Induction  problem  was  not  a  formal  part  of  the 
study,  significant  time  was  spent  on  It  because  of  Its  Importance  to  the 
problem).  It  Is  also  clear  that  the  magnetospherlc  currents  make  a  large 
contribution  to  the  day  to  day  variability  In  the  observed  pattern.  The 
fraction  of  the  observed  pattern  caused  by  the  magnetospherlc  currents  Is 
still  not  clearly  defined  for  two  reasons.  First,  there  Is  some  confusion  as 
to  what  the  magnitude  of  the  observed  pattern  really  Is;  and  second,  only 
average  values  for  the  magnetospherlc  currents  were  used  even  though  It  Is 
clear  that  for  some  quiet  days  where  the  observed  pattern  was  determined,  the 
magnetopause  current  system  Is  stronger  than  that  used  In  this  study. 

MAJOR  ACCOMPLISHMENTS  ON  THIS  CONTRACT 

The  major  accomplishments  of  the  work  completed  under  this  contract  are  listed. 

o  It  was  determined  by  direct  Integration  over  the  magnetospherlc  current 
systems  that  they  make  a  contribution  of  approximately  12  nT  to  the  day 
to  night  variation  In  the  earth's  surface  magnetic  field.  This  Is  a 
minimum  contribution  for  exceptionally  quiet  times  when  all  of  the 
magnetospherlc  current  systems  are  at  their  minimum  strength.  These 
values  are  significantly  different  from  those  given  by  using 
magnetospherlc  magnetic  field  models  directly,  since  the  models  were 
used  to  represent  the  magnetic  field  throughout  the  magnetosphere  and 
not  accurately  at  the  earth's  surface.  He  conclude  since  12  nT  Is  an 
Important  fraction  of  the  total  observed  day  to  night  variation  In  the 
earth's  surface  magnetic  field  at  mid-latitudes  that  the  magnetospherlc 

currents  produce  at  least  one-third  of  the  observed  variation  In  S  . 

Q 

This  finding  Is  Important  for  Ionospheric  physicists  working  on  the 
Ionospheric  dynamo  and  the  global  Ionospheric  wind  systems,  since  the 
winds  are  not  required  to  produce  all  of  the  observed  Sq  pattern. 
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o  It  Is  well  known  that  currents  external  to  the  earth's  surface  Induce 
currents  that  flow  In  the  earth's  crust  which  act  to  Increase  the 
horizontal  component  of  the  magnetic  field  produced  by  the  primary 
current  at  the  earth's  surface.  However,  It  Is  difficult  to  represent 
these  Induction  currents  accurately,  thus,  typically,  the  earth's 
electrical  conductivity  Is  usually  represented  as  being  a  zero  or 
Infinite  value.  We  have  developed  the  math  for  representing  the  earth 
as  a  finite  electrical  conductor.  It  was  thus  possible  to  represent 
the  total  magnetic  field  produced  at  the  earth's  surfce  by  the  ring, 
tall  and  magnetopause  current  systems.  Including  the  currents  that  they 
Induced  In  the  earth's  crust  and  their  contribution  to  the  earth's 
surface  magnetic  field. 

o  In  order  to  study  the  day  to  day  variability  In  the  pattern.  It 
was  necessary  to  develop  first  a  data  base  which  Included  observations 
of  the  solar  wind  over  a  long  period  of  time  and,  for  the  same  time 
period,  mean  hourly  values  In  the  earth's  surface  magnetic  field  at 
many  mid-latitude  observatories.  This  data  was  assembled  with  the  help 
of  NASA  Goddard  and  the  World  Data  Center  In  Boulder. 

o  Preliminary  computer  work  was  done  on  the  quantitative  determination  of 
the  signature  of  the  Blrkeland  current  system  In  Its  "Region  2"  closure 
In  the  night  side  of  the  magnetosphere.  No  model  has  been  presented  or 
published  on  this  portion  of  the  magnetospherlc  current  systems  since 
there  Is  still  some  controversy  over  the  details  of  Its  topology. 

CONTINUING  WORK 

Work  on  this  contract  has  led  to  several  conclusions  which  have  been  listed 
above  In  the  subject  on  Major  Accomplishments.  It  has  also  led  to  several 
additional  questions,  which  we  will  continue  to  explore.  Some  of  these  are 
listed  below. 

o  Work  by  several  Investigators  (especially  the  recent  work  of  Matsushita 
and  Nlshlda)  has  Illustrated  the  Importance  of  the  field  aligned 
(Blrkeland)  currents  for  variations  In  the  earth's  surface  magnetic 
field  In  the  high  mid-latitude  region.  To  do  this  work  accurately. 
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however,  we  must  first  better  understand  the  entire  complex  Blrkeland 
current  system,  Including  both  the  Region  1  and  Region  2  currents  and 
those  currents  that  flow  Into  the  auroral  Ionosphere  on  the  polar  side 
of  the  day  side  cusps.  (It  Is  also  of  Interest  to  us  In  our 
magnetospherlc  magnetic  field  modeling  work  to  better  understand  the 
closure  of  the  Blrkeland  currents  In  the  magnetosphere.  This  Is 
especially  true  of  the  Region  2  currents  which  carry  from  1  to  2 
million  amperes  over  a  relatively  confined  region  of  the  magnetosphere 
and  thus,  we  expect,  produce  a  very  characteristic  signature  In  the 
night  side  of  the  magnetosphere). 

o  We  expect  to  do  a  better  Job  on  finding  the  dependence  of  th  surface 

magnetic  field  on  solar  wind  dynamic  pressure.  Essentially,  the 

problem  to  date  Is  that  there  Is  a  lot  of  noise  In  the  data.  This 

noise  Is  typically  produced  by  the  ring  and  tall  currents  and, 

admittedly,  also  has  a  contribution  from  the  Ionosphere.  However,  we 

remain  convinced  that  the  signal  Is  there,  since  It  Is  so  large  (almost 

100  gammas  over  the  range  of  known  magnetopause  standoff  locations). 

We  have  talked  with  John  Sampson  (University  of  Alberta)  and  John  Olson 

(University  of  Alaska),  who  have  developed  filtering  techniques  which 

have  been  very  helpful  In  digging  mlcropulsatlon  signals  out  of  a  large 

noise  background.  We  expect  their  techniques  may  be  helpful  In 

developing  unambiguously  the  relation  between  magnetic  surface  field 

variability  and  dynamic  solar  wind  pressure.  If  such  a  relationship 

can  be  developed,  we  hope  to  use  It  to  determine  the  S_  baseline  and 

Q 

with  It,  an  accurate  value  for  the  earth's  main  magnetic  field. 

o  The  largest  uncertainty  In  the  work  performed  on  this  contract  relates 
to  the  earthward  extent  of  the  ring  current.  We  know  that  It  can 
change  dramatically  without  changing  the  overall  ring  current  signature 
at  synchronous  orbit  (where  most  measurements  of  the  ring  current  are 
made).  However,  the  contribution  from  the  ring  current  near  the  earth 
Is  expected  to  have  an  appreciable  signature  at  mid-latitudes  even 
during  quiet  times.  It  Is  Important  to  determine  this  contribution  to 
the  earth's  surface  magnetic  field  In  order  to  continue  the  study  of 
the  contribution  of  magnetospherlc  currents  to  S^,  and  also  for  the 
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accurate  determination  of  the  earth’s  main  magnetic  field.  We  have  had 
several  discussions  with  Suglura  and  Langel  of  NASA  Goddard  on  this 
problem. 
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APPENDIX 


THE  DETERMINATION  OF  MAGNETIC  FIELDS  PRODUCED 
AT  THE  SURFACE  OF  A  ROTATING  SPHERE  WITH 
FINITE  ELECTRICAL  CONDUCTIVITY 


The  mathematical  formalism  for  determining  the  magnetic 
field  produced  by  currents  Induced  In  the  earth's  crust  Is 
developed.  These  currents  are  formed  as  the  earth  responds 
to  external  magnetic  fields  as  It  rotates  under  them. 

Note  that  this  appendix  was  printed  with  the  TEX 
typesetting  program. 
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This  attachment  will  develop  solutions  for  the  magnetic  field  induced  interior 
and  exterior  to  a  cylinder  or  sphere  having  finite  conductivity.  The  fields  are  induced 
because  the  motion  of  the  object  in  the  presence  of  a  static  external  magnetic  field 
causes  an  electric  field  to  be  produced.  This  electric  field  will  cause  currents  to 
flow  because  the  object  has  finite  conductivity,  and  the  currents  cause  the  induced 
magnetic  field.  We  begin  with  Maxwell’s  equations  in  the  inertial  reference  frame: 


<9B 

VXE=-- 

'XI=,+i 


V  x  J  =  - 
V  -  B  =  0 


Op 

ot 


V  .  D  =  p. 


(1) 


There  are  also  constitutive  relations  that  hold  in  the  reference  frame  in  which  the 
materia]  is  motionless  (  this  reference  frame  is  denoted  by  an  *): 


B*  =  //H* 

D*  =  fE*  (2) 

J*  =  <xE\ 

The  constants  /i,  t,  and  a  are  the  material  permeability,  permittivity,  and  conduc¬ 
tivity  respectively. 

The  field  quantities  as  measured  in  the  moving  frame  -  in  which  the  material 
is  motionless  -  are  related  to  the  fields  measured  in  the  inertial  reference  frame  by 
Minkowski’s  equations  for  moving  media.  Letting  v  be  the  local  material  velocity, 
and  assuming  that  v/c  <  1  where  c  is  the  speed  of  light,  the  relations  are: 

E 
B 
H 
J 
P 
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From  equations  (2)  and  (3),  we  can  derive  two  more  useful  relations: 

e  =  e;  -  x  b 

H  =  H‘  +  rv  X  E  ' 

where  we  also  use  ep  =■  1/c2  and  neglect  terms  0(v2/c2)  as  was  done  in  equations 

(3). 

ROTATING  CYLINDER  OR  SPHERE  EQUATIONS 

We  assume  that  the  material  object  in  the  static  external  magnetic  field  is  a 
cylinder  rotating  about  its  symmetry  axis,  or  a  sphere  rotating  about  a  given  axis 
with  constant  angular  velocity  u.  In  both  cases,  the  velocity  vector  is  given  by 
v  =  uiz  x  r.  The  time  dependences  in  Maxwell’s  equations  will  eventually  damp 
out  leaving: 

V  x  E  =  0 

V  X  H  =  ctE*  +  pv 

V  •  J  =  0  (5) 

V  •  B  =  0 

V  .  D  =  p. 

Using  the  relations  in  (3),  (5),  and  V  •  v  =  — wf*  •  V  x  r  =  0,  we  can  show 
that  there  is  no  free  charge  in  the  sphere  or  cylinder.  We  derive  the  equation  p  = 
^dp/Q<p  with  the  solution  p  =  .  To  have  continuity  of  p  ,  we  must  set  po 

to  zero,  so  there  can  be  no  free  charge  in  the  body. 

The  H  field  in  the  object  can  be  shown  to  be  approximately  divergence-free  as 
follows.  From  the  V  B  =  0  equation,  we  can  show  that  V  -H—  eu)Es  =  0.  Since  u  !=» 
7.3  X  10~5  radian/sec  for  one  revolution  per  day,  and  e  «  8.9  X  10~12farad/meter, 
we  have: 

V'H«  (6.5  x  10—le amp/volt-meter)E, 

which  approximates  V  •  H  =  0  for  reasonable  sized  E  fields. 

The  V  X  H  equation  in  (5)  becomes 

V  X  H  =  irE  +  <rn(wlM  X  r)  X  H 

after  using  the  E  relation  in  (4)  and  p  —  0.  Talcing  the  curl  of  this  equation,  and 
using  the  relations  derived  above,  we  get: 

V  X  V  X  H  =  ap{u\M  x  H  —  (v  •  V)H) 

=  -‘7/‘w(^i|r-L  +  ^ + jp1’)  in  the  c?Hnder 
=  in  the  SI>bere- 


(«) 


It  is  equation  (6)  that  will  be  solved  to  determine  the  induced  fields  in  the  cylinder 
and  sphere. 


SOLUTION  FOR  ROTATING  CYLINDER 

The  solution  for  this  problem  is  found  in  a  manner  similiar  to  determining  the 
natural  modes  of  the  cylinder’s  oscillation  in  time  varying  fields.  Letting  the  H  field 
be  decomposed  as: 

H=  Y,  H„ein»  (?) 

n~  —  oo 
n 

we  find 

V  x  V  x  (Hne‘n*)  =  -inonu( or, 

V  X  V  x  M„  =  and  (8) 

V  X  V  x  Nn  =  k2nKn, 

where  k\  =  —  tnapw,  and  M„  and  Nn  are  two  independent  vector  solutions  for 
mode  n.  Stratton  (1941)  gives  the  solution  for  this  problem  by  Hertzian  vectors,  and 
by  another  method  (p  392)  which  is  also  applicable  to  the  rotating  sphere  problem. 
We  will  follow  the  latter  method. 

Briefly,  we  let  £n  be  the  solution  of 


V2f„  +  =  0. 

Then,  two  independent  vector  solutions  to  (8)  are 


(») 


Mn  —  V  X  i*£n 

Nn  =  —V  x  Mn. 


and 


(10) 


Note  that  both  M„  and  Nn  are  divergenceless.  A  solution  of  (9)  having  no  z 
dependence  is 

£»  ==:  t —Jn[knr±)etn*, 


in 


(11) 


where  J„  are  Bessel  functions  of  the  first  kind.  Therefore,  the  two  solutions  are 
M„  =  (  i-  MKr X)i,x  -  ^  J'„(knr±) and 


k 2 

N„  = 


(12) 


3C 


The  natural  modes  for  the  rotating  cylinder  combine  to  give 

OO 

H=  X)  KM„  +  6„N„)  (13) 

n  =  —  oo  ^ 

r»5^0 

where  an  and  bn  are  determined  by  boundary  conditions. 


SOLUTION  FOR  ROTATING  SPHERE 

The  solution  for  the  rotating  sphere  parallels  that  for  the  rotating  cylinder.  We 
let  H  be  decomposed  as  in  (7)  and  must  solve  equation  (8)  in  spherical  geometry. 
Stratton  has  done  this  (p  414)  by  showing  that  if  £n  is  a  solution  to  (9),  then  the 
two  independent  vector  solutions  of  (8)  are 


M„  =  Vx  i rr(n  and 

N„  =  x  M„. 


(14) 


The  solution  to  (9)  in  the  spherical  case  is 

f»=  E  C"Lfm+i(knr)P"JcoSey"*  or 
m=n  y/r 

C  =  ^m+i(*„r)P”(cosey"* 

y/r 


(15) 


where  P are  associated  Legendre  functions.  The  vector  solutions  corresponding 
to  are  and  NjJ,,  and  they  are  given  by 


sin0  d<t> 


l  Hnm. 

U - 


dO 


Nn  ^  +  l)  cW  1  9  rC,  ,  _ 1_ 

m  knr  m  r  knr  drdO  knrs\u.O  drd<p 
The  natural  modes  combine  in  the  spherical  case  to  give 


4- 


oo 


m 


h=E  E  («x+w 


m= 


11  nM0m 


(16) 


(IT) 


BOUNDARY  CONDITIONS 

The  boundary  conditions  for  these  problems  are  similiar  to  those  for  scatter¬ 
ing.  We  have  an  external  forcing  field,  a  field  transmitted  into  the  object,  and  a 
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“scattered"  or  reflected  field  externa!  to  the  object.  The  total  external  field  is  the 
forcing  plus  reflected  fields.  If  we  look  to  these  problems  as  representing  the 
variation  of  the  earth’s  magnetic  field,  it  is  the  total  external  field  that  contains  the 
S9  variation.  The  boundary  conditions  to  be  satisfied  at  the  object  surface  are:  1) 
The  tangential  component  of  the  total  external  H  must  match  this  component  of 
the  transmitted  H,  2)  The  normal  component  of  the  total  external  B  must  match 
the  normal  component  of  the  transmitted  B. 

We  have  not  talked  about  the  transmitted  B  yet  but  it  is  easy  to  show  using 
(4)  that  to  0(t/2/c2)  we  have  B  =  (iE  inside  the  object.  In  the  vaccuum  outside, 
the  relation  B  =  poH  holds  exactly. 

In  the  sections  that  follow,  we  assume  that  the  forcing  magnetic  fields  are 
caused  by  stationary  currents  that  are  bounded  away  from  the  cylinder  or  sphere. 
This  assumption  simplifies  the  mathematics  of  the  problem  because  the  forcing 
field  can  be  expanded  into  spherical  harmonics,  so  no  coupling  between  the  natural 
modes  (i.e.  of  coefficients  an  and  bn  for  different  n)  will  occur. 


* 

\ 

L 


DETERMINING  THE  COEFFICIENTS  FOR  THE  ROTATING  CYLINDER 
The  external  forcing  field  is  assumed  to  have  the  form 


(oo 

£  fUrJa^V 

n==— oo 

0 

where  a  is  the  cylinder  radius  and  the  complex  constants  dn  are  known.  The 
reflected  H  field  will  have  the  form 


B,  =  -av(  £  (19) 

\n== — op  / 

n^O 

The  tangential  boundary  condition  at  r_L  =  a  using  (12)  is 

~an{jn  =  ~dnin  ~~  gnin  (20) 

and  bn  is  zero  since  we  assume  no  forcing  in  the  z  direction.  The  radial  boundary 
condition  is 

a”(a  ' »(**“))  =  —rf^°lnl  +  0»/*o!v‘i  (21) 

The  transmitted  and  reflected  fields  can  be  found  by  solving  (20)  and  (21)  simul* 
taneously  for  an  and  gn. 
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DETERMINING  THE  COEFFICIENTS  FOR  THE  ROTATING  SPHERE 
In  this  case,  the  external  forcing  'will  have  the  form 


Hy  =  -«v(  E  E  C(’-/«)lm|C(«>^K”*) 

Vm_,  »=-m  J 

/  oo  m  v 

=  -«v(  E  E  cx") 

Sn=l  m  ' 


nj40 

with  known  complex  constants  d”  The  reflected  H  field  will  be 


•m- 
oo  m 


(oo  m  \ 

E  E  ?>A),m|+,P”(«sSK"M 

m= 1  n —  m  ' 

nj*0 

(oo  m  \ 

X)  X)  ffmC} 

m=l  B — ' 


n  j^O 

The  tangential  boundary  conditions  are 

_ 

nVaf7,sin  #  d<p 

0 


oo  m  /  t 

E  E 


+»: 


1  &rCm 


kna  drdO 


0d>" 
.f  _Z!” 

Tym  <9# 


+c^)=° 


for  the  0  direction,  and 

oo  m 


y  y  — * 

MS  msi  + 

n^O 


sin  9  3r0<fr 


sin#  dtp 


-i-  (tb  __l  aX 

l  »m_!_ /J  a  ±  rum_! 


m 


'sin#  dtp 


)  =  0 


(22) 


(23) 


(24) 


(25) 


for  the  tp  direction,  where  the  functions  are  evaluated  at  r  =  a.  Taking  (3(24 )Jd<p)— 
(i 9  sin  0  (25 )/9$),  we  get 


Y '  a*  (  1  I  AfcJn  _  ft 

~in2±m  \sin#  dtp 2  +3#^  9  90  V  °' 

fn  a  A 

AfO 
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so  0^,  =  0.  With  a JJ,  =  0,  both  (24)  and  (25)  reduce  to 


c— — 

kna 


dr 


+  CK  +  CXm  —  °' 


(26) 


The  radial  boundary  condition  completes  the  set  needed  to  determine  the  coefficients. 
Matching  the  normal  B  components,  we  get 


Ln(fitn(m  +  1) 

H  ' 


&)  = -&(<vo^)  -  Ow 


Sxl 


m 


Dr 


(27) 


Expressing  the  conditions  (26)  and  (27)  more  explicitly,  we  have 

c(  ~~  •'".+#»«)  +  =  -flC  <28) 

\2hnJa  J 

_  „•  „o(m  +  1)  =  (29) 

A[}  fl y  (Z 

and  the  coefficients  and  are  found  by  solving  (28)  and  (29)  simultaneously. 
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